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Abstract 



In order to obtain fundamental information about the Standard Model of particle physics 
from experimental measurements of exclusive hadronic two-body 5-decays we have to be 
able to quantify the non-perturbative QCD effects. Although approaches based on the fac- 
torization of mass singularities into hadronic distribution amplitudes and form factors pro- 
vide a rigorous theoretical framework for the evaluation of these effects in the heavy quark 
limit, it is not possible to calculate the 0(Aqcd/'"^6) corrections in a model-independent 
way, because of the presence of non-factorizing long-distance contributions. It has been 
argued that Sudakov effects suppress these contributions and render the corresponding 
corrections perturbatively calculable in terms of the distribution amplitudes. In this paper 
we examine this claim for the simple and related example of semileptonic 5 — ^ tt decays 
(which have similar long-distance contributions) and conclude that it is not justified. The 
uncertainties in our knowledge of the mesons' distribution amplitudes imply that the cal- 
culations of the form factors are not sufficiently precise to be useful phenomenologically. 
Moreover, it appears that a significant fraction of the contribution comes from the non- 
perturbative region of large impact parameters, and is therefore uncalculable. We also raise 
a number of theoretical issues in the derivation of the underlying formalism. Our conclu- 
sion is therefore a disappointing one. For 5-decays it is not possible to invoke Sudakov 
effects to calculate amplitudes for decays which have long-distance divergences (end-point 
singularities) in the standard hard-scattering approach. 
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1 Introduction 



The study of 5-decays is central to the development of our understanding of the standard 
model of particle physics and to the determination of its properties, parameters and limi- 
tations. The BaBar |l[ and Belle 5-factories as well as other experiments will continue 
to provide a large amount of experimental data on 5-decays in general and on two-body 
hadronic decays in particular. A measurement of the mixing-induced CP-asymmetry in 
the golden mode B — > J /'^ Kg (and related decays) for which there is a single weak phase, 
allows one to determine sin (2/5) without any hadronic uncertainties (where (3 is one of the 
angles of the unitarity triangle). For other decay modes however, our inability to quantify 
the strong interaction effects prevents us from being able to determine fundamental infor- 
mation about the standard model from the measured branching ratios and asymmetries 
with the desired precision. 

Recently it was discovered that, in the nib oo limit (where rrib is the mass of the 
6-quark), the mass-singularities in two body hadronic -B-decays factorize, so that the cor- 
responding amplitudes can be written as convolutions of universal non-perturbative quan- 
tities (the light-cone distribution amplitudes of the mesons and semi-leptonic form factors) 
and perturbatively calculable hard-scattering kernels . This provides a theoretically rig- 
orous framework for the evaluation of decays B M1M2 (where Mi and M2 are two 
mesons, at least one of which is light) in the heavy quark limit. 

The physical mass of the 6-quark is only about 4.5 GeV however, and so 0(Aqcd/"^6) 
corrections are significant (we refer to contributions which are suppressed by powers oil/rrib 
as power corrections). Indeed there are "chirally enhanced" power corrections for which the 
scale is not Aqcd but the much larger m\/{ms + md). In other cases, contributions which 
are suppressed by powers of l/nib may be enhanced by CKM factors or for other reasons. 
We would therefore very much like to be able to compute the power corrections reliably. In 
general, mass singularities do not factorize for the power corrections, and one is therefore 
reduced to using model or phenomenological estimates or exploiting flavour symmetries, 
hence losing precision and predictive power. There is a school of thought, however, which 
claims that Sudakov effects regulate these mass-singularities in such a way that they are 
calculable in perturbation theory [Q. Indeed, calculations oi B ^ M1M2 decay amplitudes 
which include power corrections and which use Sudakov effects to suppress the long distance 
effects are being presented This approach is frequently referred to as the pQCD 

formalism. 

The possibility that the power corrections may also be calculable in perturbation the- 
ory is an exciting one of course, although one may have doubts whether Sudakov effects 
are effective at scales as low as m?,. In this paper we reexamine the calculations of the 
amplitudes for a related but simpler process, B ^ ir semileptonic decays. In each order of 
perturbation theory this process is singular at long-distances. It is for this reason that the 
values of the form factors are used as a non-perturbative input in the factorization formu- 
lae 10 . On the other hand, in the pQCD approach, the Sudakov logarithms are resummed 
and the semileptonic form factors are claimed to be calculable |]^, 0. In this context we 
address two important questions: 
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1. Can Sudakov effects be used to evaluate the form-factors in perturbation theory in 
principle'! In particular we examine whether all (or almost all) of the contribution 
to the amplitudes comes from the perturbative region. 

2. Are the form-factors calculable in practicel Specifically we ask whether one can 
obtain phenomenologically useful results, given the uncertainties in the mesons' wave 
functions. 

In addition we examine some of the theoretical steps in the derivation of the pQCD formula 
for the semileptonic amplitudes in terms of the Sudakov form factors. 

The results of our study lead to disappointing answers to both questions. When we 
apply the formulae used in the pQCD approach we find that our ignorance of the mesons' 
wave functions leads to very large uncertainties in the predictions for the form factors. 
Moreover a significant fraction of the result comes from regions of phase space where 
perturbation theory (and hence the formulae which are being used) do not apply. We 
should stress that our numerical results are not in disagreement with those in refs. 0-0 
and particularly in ref. |T^. The conclusions which we draw from our calculations however, 



are profoundly different from those of these authors. We believe that our calculations show 
that the pQCD approach cannot be used to make reliable predictions for the form factors, 
and hence also in the evaluation of the B — ^ M1M2 decay amplitudes and in the subsequent 
determination of the parameters of the unitarity triangle and other studies of CP-violation. 
This is the main conclusion of our study and our motivation in writing this paper is to try 
to open a debate of this very important issue. 

We also have some reservations about the theoretical foundations of the pQCD formula 
for semileptonic decays of 5-mesons. In particular we do not understand the framework for 
the derivation of the Sudakov factor for the B-meson. We expect the Sudakov suppression 
in the 5-meson's distribution amplitude to be weak, and soft configurations with large 
transverse separations to contribute significantly to the B —>■ rriui form factors. We will 
explain our concerns in some detail in sections ^ and |^. 

The plan of the remainder of this paper is as follows. In section |^ we recall why the 
B n form factors cannot be computed using the standard factorization approach. In 
section ^ we present an outline of the framework used in the pQCD approach. The use of 
Sudakov suppression of long distance contributions has been introduced into hard exclusive 
processes involving light hadrons in order to improve the precision of the calculations, and 
is now being used also in 5-physics. In the following three sections we investigate the 
ingredients of the pQCD formalism for semileptonic 5-decays in detail and highlight our 
concerns. Section ^ contains a discussion of the distribution amplitude of the pion and 
section |^ is devoted to a critical examination of the derivation of the Sudakov factor for the 
distribution amplitude(s) of the i?-meson. We also investigate whether the approximation 
of using a single distribution amplitude for the 5-meson is a valid one. In section ^, 
the pQCD expressions for the 5 — > vr form factors are shown to be very sensitive to the 
shape of the mesons' distribution amplitudes, and to contain significant uncalculable long- 
distance contributions. Section |^ contains a summary of our principal conclusions. In the 
appendix we introduce a model distribution amplitude for the 5-meson, which satisfies 
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Figure 1: Diagrammatic representation of the i? — > tt semileptonic decay in the factorization approach. 
The thick hne represents a 5-quark. 



the constraints due to the equations of motion in longitudinal and transverse momentum 
space. 



2 B ^ 7T form factors in the factorization approach 

In this section we consider the semi-leptonic decay B{p) — > 7r(p')£z/£ in the standard fac- 
torization approach, introduced by Brodsky and Lepage and by Efremov and Radyushkin 
(BLER) [|lT], |12|. We introduce the main ingredients of this approach, and in particular the 
light-cone distribution amplitudes which describe the parton content of the mesons relevant 
for this process. We then recall why long-distance singularities prevent the computation 
of -B — i> vr form factors in the standard approach [l^ . 

We start by defining the kinematics and our notation. The amplitude for this decay 
can be written in the form: 

Aip,p') = ^ K6 - ( vr(p') I u^,b I Bip) ) , (1) 

where / and ui represent the wave functions of the charged lepton and neutrino respectively 
and u and b are quark fields. The non-perturbative QCD effects are contained in the 
hadronic matrix element ( tt{p') \ u'~f^b \ B{p) ), which can be parametrized in terms of two 
invariant form factors. We choose to use a conventional definition of the form factors (in 
a hehcity basis): 



; n{p') I uj,b I B{p) ) = F^{q'){p + p'), + (M| - Ml)^^^^^^-^^q, , (2) 



where = p^—p'^. At = we have -fo(O) = -F+(0). The notation is illustrated in Fig. |I]. 

We work in the rest frame of the i?-meson, and neglect the mass of the pion. We 
introduce light-cone coordinates, writing k = (A;+/-\/2, A;_/-\/2, with k± = (fco ± ^3), 
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and with the metric for the transverse components chosen such that {k±)^{l±)fj^ = —k± ■l±. 
It will also be convenient to introduce two light-like vectors n+ = (-\/2, 0,0_l) and ?2_ = 
(0,72, Ox). 

Using light-cone coordinates, p = (Mb/ V2, Mb/ V^,0±) and we chose the momentum 
of the pion to be in the — direction, p' = (0, tjMb / v^, 0^)- In terms of r], the energy of the 
pion is given by p ' ° = r^M^ /2 and the invariant mass of the lepton pair by m|j^ = (1— ?7)M^. 
The physical range for rj is therefore given by < r/ < 1. 

In the standard hard-scattering approach, as represented in Fig. |l|, the form factors are 
obtained from convolutions of the mesons' distribution amplitudes (0,r and 0^) and the 
perturbative kernel T. In the following subsection we discuss the distribution amplitudes 



before returning to the convolution itself in section 2.2 



2.1 Meson distribution amplitudes 



In order to evaluate the decay amplitude, we have to perform the convolution illustrated 
in fig. |I|. In this subsection we discuss one of the principal ingredients of this convolution, 
namely the distribution amplitudes (or DAs) of the incoming and outgoing mesons. We 
will consider only the leading-twist amplitudes here. 



The (leading-twist) distribution amplitude of the pion is defined as |jTT 



{■^ip')\qaiy)q'(iiz)\o) = i^{^'j5)(3a dxe 



i{xp' -y+xp' -z) 



(3) 



where {y — z)'^ = 0, x = 1 — x, and we choose a normalization such that — 131 MeV. 
A path-ordered exponential is implicit on the left-hand side, in order to maintain gauge 
invariance. The parameter /i is a renormalization scale of the light-cone operators on the 
left-hand side. 

The distribution amplitude is symmetric with respect to the interchange x ^ x, and is 
normalized such that /q^ rfx 07r(x; /i) = 1. As the renormalization scale /x tends to infinity, 
0^ tends to the asymptotic limit 0^^(x) = 6x{l—x). For finite values of the renormalization 
scale, can be expanded on the basis of Gegenbauer polynomials: 



07r(x; fi) = 6x{l — x) 



l + E«2p(/i)Cg/')(2x-l) 

p=l 



(4) 



Only even polynomials contribute because of the symmetry x ^ x. The lowest Gegenbauer 



polynomials contributing to the sum are C2 



(3/2), 



U 



3/2 ■ (5^2 - 1) and a 



(3/2), 



U 



15/t 



{21u^ — lAv? + 1). The moments are multiplicatively renormalized, with the anomalous 



dimensions increasing as n increases [14| (so that the corresponding a2p decreases more 
quickly with increasing p a.s fi 00). 

For the distribution amplitudes of the i?-meson, we will use the notation and formalism 
developed in refs. and [0 which we now briefly summarize Q. In the heavy-quark limit. 



^We note the reservations concerning the use of such distribution amplitudes for heavy mesons at next- 
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the most general decomposition of the light-cone matrix element for the 5-meson is |]T5|, p!6[ : 



OMz)b^{0)\Bip)) 



4 



t 



75 



(5) 



«/3 



where z'^ = 0, v = p/m, t = v ■ z and a path-ordered exponential is implicit on the 
left-hand side again. We are working only to leading twist order and therefore have set 
Mb = nih = m. Recalling that ( 0|g(0)7^756(0)|i?(p)) = ifsP^, the choice of prefactor in 
eq. implies that 

0^(t = O)=0^(t = O) = l. (6) 

Thus for the 5- meson we have two distribution amplitudes, 0^, to consider. 

In order to evaluate the amplitudes for B iriue decays we perform a convolution of the 
distribution amplitudes of the 5-meson with the hard-scattering amplitude, as illustrated 
in Fig. |1|. Let / denote the momentum of the light quark in the 5-meson and recall that we 
have chosen the momentum of the outgoing pion to be in the — direction. The 0{as) hard- 
scattering amplitude T is then independent of We are thus allowed to set z^. = in the 
decomposition in eq. (^). Following ref. |jl6[ it is possible to write the convolution in the 
form dl+ {P^T(l)}, with / set to {u /\/2,0,0±) once the momentum-space projection 
operator has been applied. The projection operator is equal to: 



pB 



(7) 



a/3 



with A(co') = l'^ di — 0+(^) )• The distribution amplitudes on the right-hand side of 

eq. (|^ are written in momentum space. 



<Pl{uj) = -j_Jte^-'m)- (8) 

Since the components of the light quark in the S-meson are expected to remain of or- 
der Aqcd or less, the distribution amplitudes will be suppressed for ui larger than Aqcd- 
Note also that eq. (|^) implies that the distribution amplitudes in momentum space are 



normalized such that / duo < 



J duj (j) 



1. 



For a precise calculation of the form factors it is clearly important to have as much in- 
formation about the mesons' distribution amplitudes as possible. Of particular importance 
is the behaviour of the distribution amplitudes at the end-point /+ — > 0. To gain some in- 
sight into this behaviour, we follow ref. |TB[ and note that the leading-twist decomposition 
in eq. leads to expressions for the distribution amplitudes in terms of matrix elements 
of the pseudoscalar density and axial current: 



{0\q{z)j,b{0)mp)) 
(0 I g(^)7M75&(0) \Bip)) 



-^IbMb 
Ib 



^1, 



and 



(9) 

(10) 



to-leading order in the perturbative expansion expressed in refs. p7| , p^ . At this order, which is beyond 
the scope of our study, it is hkely that the formahsm will have to be modified, at least to include the 
dependence on the transverse momentum explicitly. 
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where 



and 



are given in terms of the distribution amphtudes (f)^ by 



lB 



^1 



+ 



kB 



^^1 



and 



'.B 



t 



'111 



For a hght pseudoscalar meson, the end-point behaviour of the eigenfunctions of the evo- 
lution equations, together with sum- rule inspired arguments, suggest that (j)^-^^{uj) ~ u and 
4>p{u}) ~ 1 for u; — s> [|12|. If this behaviour is also valid for the 5- meson, we might expect 
that = and = 0(1) as /+ ^ [1^. However, in spite of this argu- 

ment, when studying the validity of perturbative estimates of the form factors in terms of 
the distribution amplitudes it must be remembered that such an end-point behaviour is a 
conjecture. 

Another theoretically motivated constraint is a consequence of the equations of motion. 
Since we neglect the 3-particle distribution amplitudes (and contributions from higher Fock 
states), the equations of motion lead to the constraint |T6[: 



dL 



(12) 



Eq. (0) is useful in constraining models for the distribution amplitudes, and supports the 
conjectured behaviour of the distribution amplitudes as /+ — > 0. We can indeed derive 
the expected behaviour, that 4>^{l+) = 0(1) as /+ — >• and that 0+(O) = 0, if we assume 
that the equation of motion ([T2|) is fulfilled and that A^^ = J^dl+ is finite and 

nonvanishing 0. We can use the equation of motion to see that A^^ = 0:^(0) [|^, which 
leads to 0:^(/+) = 0(1) for small /+. Eq. ([T2|) implies then that 0:^(/+) vanishes at the 
end-point. 

We end this subsection by introducing two models for the distribution amplitudes which 
satisfy the above constraints, and which will be useful in our investigation of the validity 
and reliability of the frameworks used to evaluate the form factors. The first model was 
proposed in ref. on the basis of a QCD sum rule analysis: 



kB 



A 



= -exp 



A 



(13) 



A is a parameter of order Aqcd- Another model satisfying the above constraints is discussed 



in Appendix A.l 



2 / 



2A2 



ttA' 



exp 



2A2 



(14) 



In this model, as discussed in the appendix, A is related to the distribution of transverse 
momenta. 
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Figure 2: Diagrammatic representation of the lowest-order hard-scattering kernel for B — > iriiye. 
transverse momenta are introduced for the discussion in sections ^ and ^. 



The 



2.2 Expressions for the B ^ tt form factors 

We are now ready to discuss the semileptonic B ^ tt form-factors in the standard hard- 
scattering approach. The leading-order 0{as) contribution to the i? — > tt form factors is 
given by the two diagrams in Fig. ^ and leads to the expressions: 



;.;,;!=^^',,/,,,Mli_g+Ml, (15) 



The expression for appears in ref. [0. In deriving these expressions we have used the 
constraint from the equation of motion eq. ([12D to write A{uj) = Jq di {(t)^{€) — 0+(^) ) = 

From equations ( plSf ) and (|16D we see that the B —>■ tt form factors are very sensitive 
to the end-point behaviour of the distribution amplitudes for the pion and the B meson. 
If we suppose, for example, that (pni^) ~ x for small x and (p^il^) ~ 1 for small /+, the 
integrations over x and /+ are both divergent. It is therefore not possible to compute 
B 71 form factors in the standard BLER approach. It is for this reason that the authors 



^We thank Martin Beneke for raising this point. 
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of ref. proposed to consider these form factors as non-perturbative inputs (obtained for 
example from experiment or lattice QCD) when evaluating two-body hadronic 5-Decays. 
In the pQCD approach on the other hand, Sudakov effects are invoked to regulate the 
integrals at the end-points. We now proceed to discuss this approach. 



3 Outline of the pQCD approach 

A different approach (modified factorization or pQCD) has been introduced with the aim 
of solving the problem discussed in section ^]2] above. In refs. it is claimed that the 

introduction of transverse momenta and the subsequent resummation of large logarithms 
into Sudakov form factors regulates the end-point singularities sufficiently to enable the 
B ^ 71 form-factors to be computable. In this approach, the form-factors are expressed as 

Fo+= / dx I dl+ d'^b±d'^c±'^^{l+, c±; fi) Hq +{xj+,b±, c±,r],m; ^)'^Tr{x,b±-, fi) , (17) 
Jq Jo J 

where: 

• the meson distribution amplitudes '^■j^ and '^^ contain the resummed Sudakov effects; 

• b± and c± are the variables conjugate to the transverse momenta of the valence quarks 
and correspond to their transverse separations; 

• i^o + are the Fourier transforms of the corresponding hard-scattering kernels. 



In the following sections we will examine the ingredients in eq. (|T^ in some detail. For 
presentational purposes however, we hope that it will be instructive if we begin with a brief 
discussion of the philosophy of the pQCD approach in the simpler and much-analyzed case 
of the pion electromagnetic form factor This form factor can be computed in 

the standard approach |[Tl|, ^ , leading to a convolution of a hard-scattering kernel T with 



the DAs of the two pions: 

F^--7*(Q)= f' dudvMu;fj.)T{u,v,Q;fi)Mv.l^) , (18) 
Jo 

where Q is the momentum transfer. Although the convolution in eq. (18) is finite, it has 



been argued ||2^, ^ that substantial contributions may come from end-point configurations 
{u,v — >■ 0, 1), where one or more propagators of the hard-scattering kernel T gets close to 
the mass-shell, enhancing the soft contributions. It is likely that very large momentum 
transfers may have to be reached in order to ensure a sufficient suppression of these soft 
contributions. 

The authors of ref. ||20| proposed to exploit the ideas developed in refs. |2^, ^ in 



order to extend the domain of applicability of eq. (|T^). The essential idea is to keep the 
transverse momenta of the valence quarks writing: 

F^^^^* = du dv j d^k^ dH^ ^^(u, k^.] /i) H{u, v, k±, Q; fi) "^^{v, l±] /i) . (19) 
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Eq. ( p!9D is analogous to the Fourier transform of eq. (jlT]). The tree- level contribution to the 
hard-scattering kernel is now modified. A typical modification is l/{uvQ'^) l/luvQ"^ + 
(A;_L + thus using the transverse momenta to regulate the end-point singularities. 

Since the transverse momenta are soft, further arguments have to be invoked in order to 
demonstrate that the procedure is a consistent one. 

The modified hard-scattering kernel is now convoluted with a light-cone wave-functions 
which depend on both longitudinal and transverse components of momentum. Such wave- 
functions are defined by [^ : 

lUp'^^^ix, k^) = J ^|^e^(^^'-^--^^-"-)(0|g(^)7.75g'(0)|7r(p'))L,=o (20) 

in an axial gauge n ■ A = [n^ 0]. We emphasize that no path-ordered exponential is 
present in the matrix element of eq. (^); "^t, is a gauge- dependent quantity (indeed, as 
explained below, the gauge dependence is exploited to simplify the calculation). 

It is argued that gluon exchanges suppress configurations where the quark-antiquark 
pair has a large transverse separation. Soft and coUinear gluon exchanges between the 
valence quarks lead to large double-logarithmic effects, which need to be resummed. This 
resummation is performed in impact parameter space and so we define the Fourier conjugate 
of the wave-function, 

^n{x,b^) = f ci2fc±e^^^-^^^(x,fcx) • (21) 



Expression (|I9|) for the form factor can be re-expressed as an integral over the impact 
parameters. 

An analysis of the double logarithms from higher-orders of perturbation theory shows that 
they can be resummed and included by multiplying the wave-function in impact parameter 



space by a Sudakov factor p^. In the leading-logarithm approximation 

'^.(x, 6±) ^^^(x, Q; /i) = exp[-5^^(6^, g)]0^(x; /i = 1/6^) , (23) 

where 

S^^{b^,Q) = ^\n\b^Q) . (24) 
47r 

When subleading logarithms are included, S becomes a function of log[a;Q/AQCD] and 
log[l/(6_LAQCD)] and diverges when either of these logarithms vanish (see eq. ( ^4]) in the 
following section). This factor suppresses two regions: large transverse separations b± 
{b± ~ 0(1/Aqcd)) and small fractions of longitudinal momentum a; (x ~ 0(Aqcd/Q))- 



For the electromagnetic form factor of the pion the Sudakov factor is considered only 
when it suppresses the configurations with a large transverse separation, and it is set to 1 
outside of the region of large b±. Potential suppression for end-point configurations (small 
x), and (small) enhancements for intermediate x and b± have been neglected. 
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What is generally taken for the scale of the coupling constant in the hard-scattering 
kernel H in eq. (0) is the largest available virtuality for the hard gluon, i.e. the maximum 
of s/wvQ (longitudinal), and 1/c_l (transverse). It is then argued that the coupling 
constant remains small and the perturbative approach is consistent, unless h± and c±_ are 
both 0(1/Aqcd) and m or f — > 0. Moreover, although the coupling constant is large in 
this soft region, the corresponding contribution is suppressed by Sudakov factors of the 
distribution amplitudes. The contribution from the soft region is therefore argued to be 
small. 

Before returning to the semileptonic B —* it decay we briefly summarize the preceding 
paragraphs about the possibility that the modified factorization approach might improve 
the BLER approach. Its principal aim is to treat the end-point regions of phase space, 
in order to improve the accuracy of perturbative calculations. The singularities of the 
propagators are smoothed by considering transverse components of the momenta. The 
scale of the hard-scattering kernel is taken to be the maximal (longitudinal or transverse) 
virtuality of the exchanged gluons. The only region of large corresponds to small 
longitudinal momenta and large transverse separations. This region is expected to be 
strongly suppressed by the Sudakov exponential and it is therefore argued that a consistent 
perturbative computation is possible. 

The authors of refs. 0, |T^ claim that a similar procedure can be followed to compute 
B ^ 71 form factors within the pQCD approach. Let us first notice that this process is 
qualitatively different from the transitions vr — TC'-f* [^-[^ and 77* ir^ |2^. For the 
latter form factors, the BLER approach yields finite answers in the large limit, which 
might however be affected by significant corrections, even at large energies. The pQCD 
methods are then designed to improve the accuracy of the computation, and to investigate 
lower momentum transfers. On the other hand, we have seen that B ^ it form factors 
cannot be computed in the usual factorization approach, since singular long-distance effects 
arise. The pQCD approach would therefore not only improve the precision of the BLER 
framework, but more importantly it would regulate the divergences. In refs. [^-[0 it is 
claimed that indeed accurate and reliable results can be derived with no breakdown of the 
perturbative framework. 

The main motivation for this paper is to examine the basis for such a strong claim. In 
the following sections, we examine the various elements involved in eq. (p!7|), and emphasize 
our concerns about the validity and reliability of the pQCD approach 0. 



•^The pQCD approach that we have outhned is referred to as "fc_L-resunimation" . Recently, it has 
been proposed to combine it with a second ( "threshold" ) resummation to study i?-decays |29| . The 
concerns we express in sections ^ and || about the pQCD approach do not depend on whether the threshold 
resummation has been performed or not, and we have checked that our numerical analysis in section ^ is 
almost unchanged when both resummations are considered. We therefore focus on the fc^-resummation in 
the remainder of this article. 
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4 Sudakov effects in the distribution amplitude of the 
pion 



In this section we consider the first ingredient of eq. ([17|), the distribution amphtude of the 
pion ^^Tr, including Sudakov effects. The pion's momentum is taken to be in the — direction, 
■p' = {0,Q/\/2,0±) (where the components are in the +, — and transverse directions 
respectively), and we consider a configuration in which the valence quarks have longitudinal 
momenta xp' and xp' and a transverse separation b. The authors of ref. |^ obtain the 
following expression for the distribution amplitude "^t^ in an axial gauge [n- A = 0, n"^ ^ 0]: 

^^{x, b, Q- /x) = exp [-S{n, Q, b)] exp [-^(/x, b)] 0,(a;; /x = 1/6) + 0(a,(l/6)) , (25) 

where: 

• iS is the Sudakov factor, which contains the resummed double logarithmic contribu- 
tions, as well as some non-leading logarithms which are expected to be significant. 
It depends on the transverse separation b, on the pion's momentum p' and on the 
gauge-fixing vector n; 

• £^ is an evolution factor, relating the renormalization scale /i and 1/6; 

• is the standard BLER distribution amplitude. 

Eq. (p5|) exhibits two important features. The first one is the gauge dependence of the 
Sudakov term S. Since the convolution in eq. (|T7p yields physical (gauge-independent) 
form factors, the gauge dependence of the distribution amplitude "^.^ has to be cancelled 
by that in the hard-scattering kernel. The second issue is the presence of 0{as) corrections 
in eq. (|25|). The modified distribution amplitude can only be replaced by the standard 
one (pTT (as in eq. (^) ) for small transverse separations b. Specifically, there are relative 
0(as(l/6)) corrections to eq. (^3]). Although the Sudakov factor should suppress the full 
distribution amplitude at large b, if this suppression is weak (or negligible), the expression 
(^) can not be used reliably in this region. We investigate this numerically in section ^ 
below. 

We will return to these issues later, when we consider i? — > vr form factors. We now 
continue with the discussion of Sudakov effects in the pion's distribution amplitude and 



summarize the analysis performed in ref. |2^, and then discuss the gauge- dependence of 
the Sudakov factor. 



4.1 Derivation of the Sudakov factor for the pion 

Sudakov effects are expected to suppress configurations with large transverse separations. 
The explicit expression arises from the resummation and subsequent exponentiation of 
double logarithms. Such logarithms are related to the exchange of soft and coUinear gluons 
and double logarithmic contributions are due to an overlap of both types of divergence. 
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Axial gauges n ■ A = (where is a fixed vector in the (+, — ) plane and 7^ 0) 
are particularly suitable for use in studies of Sudakov effects. In these gauges, the overlap 
arises when a gluon is exchanged between the valence quarks of the same meson [two- 
particle reducible diagrams]. Sudakov effects can therefore be analyzed independently of 
the physical process, and are included fully in the distribution amplitude as can be seen 
in eq. (pSf). It would however be reassuring to have an explicit demonstration of the gauge 
invariance of the form-factor in which the gauge- dependence of the distribution amplitudes 
is shown to cancel that of the hard scattering kernel. 



In spite of the technicalities it is instructive to recall the derivation of eq. ( P5|) |jTO|, ^ 



22| , Pq , p7| . As mentioned above, the distribution amplitude in x and impact parameter 
space, \E'^, is defined in an axial gauge, and is gauge dependent. The only large invariant 
parameter that can be constructed in the axial gauge is z/^ = {p' ■n)'^ /\n'^\. The distribution 
amplitude is therefore a function \l/7r(x, 6, Q; /i), where the Q-dependence is linked to the 
gauge dependence, and /i is a renormalization scale. As we shall see below, the gauge 
dependence can be used to determine the Q-dependence. 

Eq. (P5|) is obtained by studying the Q-dependence of '^T,{x,b,Q; fi). This distribu- 
tion amplitude is expected to exhibit the structure '^T,{x,b,Q; fi) = exp[—S{x,b,Q)]x 
\E'^(a;, 6; /i), where \E'° is independent of Q, and S is found to be a RG-invariant function 
containing double logarithms of Q. Taking such an exponential form as an ansatz, we have 

S contains at most double logarithmic terms in Q. QdS/dQ therefore contains at most 
single logarithms of Q, which can be treated by usual renormalization-group methods. 

To obtain the pion's Sudakov factor, we have therefore to find a differential equation of 
the form, Qd^^^^/dQ = C ■ \E',r (where C contains single logarithms of Q) and to integrate 
this equation. The procedure can be summarized in the following steps: 

• \E',r depends on Q only through u"^ = {p' ■ n)^/|r2^|, so that 

^ dQ ~ p'-J dn^ • ^ > 

• The n-dependence of Feynman diagrams is contained only in the gluon propagators. 
In the axial gauge, the propagator of a gluon with momentum q is 

r,^^(,)^ZlL^^^l±^^n^jtl^, (28) 
q \ n ■ q [n ■ qY J 

and the derivative with respect to n" can be rewritten in the form: 

^ -N^'^iq) = ^(Ar^"g^ + N'^'^q'') . (29) 



dua q ■ n 
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Figure 3: Diagrammatic representation of the application of eq. (|29| ) for the derivative of a gluon prop- 
agator with respect to rt". The dashed arrow represent the injection of momentum q, using the standard 
vertex —igjpT"^ multiphed by q'' {p — n or v). 




Figure 4: Diagrammatic representation of the derivative of ^'jr with respect to logQ. The square vertex 
corresponds to —igT^ {p'-n)(q-n) P'°' ■ 



• Consider a Feynman diagram contributing to "^t^. The derivative with respect to 
logQ is equivalent to modifying one of the gluon propagators according to eqs. ( pT]) 
and (pQl). This is represented diagrammatically in Fig. ^ for the case of a gluon 
propagator attached to quark lines. We need to compute the effect of Qd/dQ on all 
the diagrams of order 0(a"), and so we have to consider all the diagrams of order 
0(a"~^), inserting in all possible ways two elements: an external gluon propagator 
carrying momentum q (wavy line in Fig. ^ and an injection of momentum q (dashed 
arrow in Fig. An integral over the momentum q has to be performed. 

• When we sum over all the possible momentum injections (represented by the dashed 
arrows in Fig. Ward identities lead to cancellations. The only surviving graphs 
correspond to a momentum injection at the end of one of the valence quark lines. 
Thus, applying Qd/dQ means inserting a gluon line, with one leg attached to a 
valence quark line with the vertex (denoted by a square in Fig. H): 

[p' ■ n)[q ■ n) 
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For the modified diagrams, the leading contribution for large Q arises when the 
inserted gluon is either soft or ultraviolet. The corresponding terms have been com- 
puted at the leading order in as in ref. yielding the differential equation: 



Q-^'^nix,b,Q;ix) 



(31) 

where /C and Q collect respectively soft and ultraviolet contributions, and have op- 
posite anomalous dimensions 7ac = — 7g- The bracketed factor in eq. ( ^I]) is therefore 
scale-independent (in ref. it is shown that this independence of the scale is also 
true at higher orders). The explicit expressions (at order 0{as)) for JC and Q in the 
MS scheme are 



g{xu/fi,g{fi)) 



and 1 26, 30 



Ik = -76 



3 TT 



+ 4 



67 
18 



4a, 

'3 TT 

4a, 
'3 TT 



7r 



log(6V'e'V4) 
2 log 



N 



6 27 



+ Oial 



(32) 
(33) 

(34) 



The differential equation eq. (plf ) can be solved directly, since the Q-dependence of 
}C{bfj,, g{fi)) + Qi^xu/ fi, g{fj,)) comes only from the i/-dependence of Q. A different 
method was followed in refs. |2^, ^ . Using RG-equations, all the Q-dependence can 
be absorbed into the scale of the coupling constant: 



}C{bix, g{ix)) + g{xiy/ IX, g{ix)) 



/C(Ci, g{C2Xu)) + g{l/C,, giC^xu)) (35) 
-2 



C2XU fjii 

^A[C,,gm)] 

Ci/b fi 



where Ci and C2 are free parameters, which should be tuned to avoid large logarithms 
in the perturbative expansion of JC and Q on the right hand-side of eq. (^). The 
function A is given by 



A[C,,g] 



lK{9)+m§-glC{Ci,g) 



Combining eqs. (0) and (|35D yields the general solution: 

*^(a;,6,(5;/i) = exp[-5(a;, 6, 1/)]^^° (x, 6; /x) , 
S{x,b,h') = s(x, 6, z/) + s(x, 6, z/) , 
1 



s{x, b, v) 



rCiXU (In 


2 log 




Jci/b p, 







AiCugifi)) 



1 



(36) 



(37) 
(38) 

(39) 



.}C{C^,g{fi))-g{—,gm) 
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We can rewrite S more conveniently using eqs. (|3^), fp^ 
expression for the strong coupling constant: 



1 



6|), and the two-loop 

/?2l0gl0g(/iVA2jCD) 



TT 



In this way we obtain 



A log(/xVA2jcD) Pf log^(/xVA|cD) 
33 - 2A^/ 153 - 19A^/ 
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/32 



24 



5(a;, 6, z/) = s{x, b,r ■ Q) + s{x, b,r ■ Q), 



(40) 
(41) 

(42) 



where r = y\X/p\ depends on the components of the gauge-fixing vector n which 

we write as n = (A/a/2, p/\/2, 0^). The gauge dependence now manifests itself as a 
dependence on r rather than u. s{u, b, R) is a function of: 



q = log 

and is given explicitly by: 

AW 



C2UR 



2A, 



QCD_ 



and b = log 



bA 



QCD 



(43) 



s{u, b, R) 



q 



2/3i 



glog U + 



f-1 



Ad) 



(q-b) 



log(26) + 1 log(2g) + 1 



4/?f 4/3: 



log 



W 27-1 
CI 



log 



9 

'q" 
.b> 



(44) 



8/3? 



log2(2g)-log2(26) 



The coefficients /?j are defined in eq. (^T]) and the A^^^^'s in eq. (|4^) are given by 

.2 in s ^Cxe^ 



4 
3 



67 



TT 

y 



10,, 8^ , 
^iV, + -/?,log 



(45) 



These coefficients come from the two-loop expression of A defined in eq. (|36|) : 



TT V TT / 



(46) 



The first three terms of J^'^^ come from the two- loop anomalous dimension of /C , 
and the last one from the partial derivative of /C with respect to g. 

Eq. (|4^) has been derived using the two-loop level expression of the coupling constant 
and of 7x:, but the one- loop expression of the functions /C and Q . It is not necessary 
to compute their two- loop expressions if C\ and C2 are tuned in eq. (39) in order to 
avoid large contributions of 0{c?^. A conventional choice is Ci = 1 and C2 = 

p^, although we are not aware of an argument that this choice is optimal (or 



even good enough for the precision claimed for 5-decay form factors p[]-|p.U|]). 



15 



Finally we need to consider the residual function in eq. (p^). Since the 

Sudakov factor suppresses the modified DA for large b, we have, in particular, to 
determine 6; fi) for small b. It is argued in ref. |27| that the factor exp{ik± ■ b) 

in the integrand of the Fourier transform in transverse space ensures that the integral 
is dominated by the region \k±\ < 1/6, and hence that 



^^{x, b;fi= 1/6) = /i = 1/6) + 0(a,(l/6)) . 



(47) 



The logarithmic corrections of 0{as{l/b)) arise because the way that the integral is 
regulated at large transverse momenta is different from that in the standard definition 
of the BLER distribution amplitude. A second type of correction is due to the 
different definition of the BLER and modified DAs. Both quantities are defined from 
the non-local, and gauge dependent matrix element (0|g(2;)7^75g'(0)|7r(p)) in a gauge 
n ■ A = 0. But (pT, is defined in the light-cone gauge [n^ = 0], whereas "^j, is defined 
in an axial gauge [n^ ^ 0] . The identification of 0,^ and \E'7r holds only in the leading 
logarithmic (LL) approximation. Thus there are next-to-leading logarithmic (NLL) 



corrections in Eq. (^T]). 

The presence of these corrections has important implications for the phenomenolog- 
ical applications of this formalism, as will be discussed in detail in section ^ In 
particular, for the calculations to be valid, the integrals have to be dominated by 
the region of small 6, so that 0^(1/6) is in the perturbative regime. The identifica- 
tion of "^^{x, 6; /i) with 07r(x, fi) is valid for /i = 1/6 and so we have to perform the 
RG-evolution of "^t^ to the scale fi = 1/6. We finally arrive at: 



"^Trix, 6; /i) = exp 



^l /i 



$^(x,6;l/6) 



(48) 



where 7^ = — as/vr + 0{a1) is the anomalous dimension of the quark's wave function 



in the axial gauges . 



Combining eqs. (pZD , (^Tf) and (^Sf ), the distribution amplitude \I'^ is the product of three 
factors: a Sudakov exponential, an evolution-related term and the standard distribution 
amplitude (p.,,. The combined expression is: 



'^^{x,b,Q;fi) 



e exp 



0^(a;;/x = l/6) + O(a,(l/6)) . 



(49) 



4.2 



Dependence of the Sudakov form factors on the gauge and 
on the integration constants. 



In this subsection we add some comments on the dependence of the Sudakov factor eq. (44) 
on Ci, C2 and n. The general expression was introduced in ref. but in the subsequent 
papers |l^, ^ , only the particular expression withCi = 1, ^2 = V2and n oc (1, —1, ) 
was considered. 
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Ci and C2 are arbitrary parameters, introduced in order to solve the differential equa- 
tion eq. (|3lD . The full expression for s{u, b, R) is therefore independent of these parameters. 
However, eq. (0) is a truncated expression, derived from the one-loop expression of /C and 
Q and the two- loop expression of 7^;. It is therefore weakly dependent on Ci and C2. A 
redefinition of these two coefficients does not affect the double logarithms (i.e. the terms 
proportional to log(g/6)) but will have an effect on the single logarithms. Ideally Ci and 
C2 should be chosen so as to minimise the effects of the unknown higher-order terms. 
This requires some insights into the important regions in x and b. The arguments of the 
logarithms in eq. (|^) are positive for the regions 

, , C*! , 2Aqcd , , ^ 2Aqcd 

and Sudakov suppression arises as b reaches Ci/Aqcd- Let us choose a specific gauge (r) 
and a momentum transfer {Q). If we consider a soft configuration for the valence quarks 
of the pion, the estimated Sudakov suppression varies with the choice of Ci and C2. Of 
course this dependence on Ci and C2 is due to the fact that we have truncated the series 
in as to obtain eq. (|4^). 

\I'^ was explicitly defined in a gauge-dependent way in eq. (^O]). As we have seen 
above, this gauge dependence can be used to determine the functional dependence of the 
distribution amplitude on Q. The gauge dependence is present in eq. (^9[) through the 
Sudakov term S. To exhibit this gauge dependence more transparently we use our freedom 
to redefine the constants of integration Ci and C2, writing Ci = C( and C2 = C'j/'", so 
that 

where the ellipsis denotes single logarithms. In eq. (0) we see explicitly the gauge depen- 
dence of the Sudakov factor. 

Of course the convolutions which yield physical form factors, such as eq. (p!9D, are gauge 
independent. The gauge dependence of the distribution amplitudes has to be compensated 
by a similar term in the hard-scattering kernel. For the case of the hard elastic scattering 
of two mesons, the authors of ref. explain clearly how the resummation of soft gluon 



interactions in the hard-scattering kernel leads to a gauge dependence cancelling the one 
in the meson distribution amplitudes. We are not aware of similar explicit demonstrations 
for the electromagnetic form factor of the pion or for the B ^ tt semileptonic decay form 
factors. Such demonstrations would add confidence that all Sudakov effects are correctly 
included in the mesons' distribution amplitudes. 



5 B-meson distribution amplitude (s) 

In section we have seen that there are two distribution amplitudes for the 5-meson 
at leading twist. In the applications of the pQCD framework to 5-decays 0-||lO|) oiie 
linear combination of the two distribution amplitudes is neglected; specifically the two 
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distribution amplitudes are set equal to each other, = = "^^^ We consider the 



validity of equating the two distribution amplitudes in section The remaining B- 
meson distribution amplitude is analyzed along similar lines to that of the pion. Consider 
a i?-meson of momentum p = {m/ \^,m/ ^/2,0±) containing a light quark of momentum 
/. The modified distribution amplitude used in pQCD analyses is: 

vI/^(/+,6,p; /i) =exp[-5'(n,p,6)]exp[-^(/i,6)]0^(Z+;/i= l/6) + 0(a,(l/6)) , (52) 

where: 

• S' is the Sudakov factor containing the resummed double logarithmic contributions 
as well as selected nonleading logarithms. It depends on the transverse separation b, 
on the momentum of the 5-meson p and on the gauge-fixing vector n. 

• £^ is an evolution factor, relating the renormalization scale fi and 1/6. 

• (f)^ = (j)^ = (f)^ is the "common" distribution amplitude of the 5-meson. 

The presence of the heavy b quark significantly modifies the analysis of Sudakov effects. 
The exponentiated double logarithms of the Sudakov factor sum the overlapping contri- 
butions of soft and coUinear gluons. We work however, in the rest frame of the B-meson 
so that the term "collinear" is not defined until the convolution for the form factor is in- 
troduced (with a hard scattering kernel which is independent of The definition of the 
distribution amplitude therefore cannot be completely separated from the process being 
studied. Moreover, all the components of the momentum of the light quark are O(Aqcd); 
so that the kinematics is very different from that of an energetic pion with a light-like 
four-momentum. In view of these questions we examine the derivation of eq. (0) in sec- 



tion 5.2 below. We start however, by considering whether setting = is theoretically 



motivated. 



5.1 Two distribution amplitudes or one? 



In refs. [g-|IO[ the phenomenological analysis is performed with a single distribution ampli- 
tude for the B- meson, setting (j)^ = (p^ = 0^ and = = \E':^. Such an identification 
simplifies considerably the calculation of the form factors and other physical quantities, 
since the momentum-space projection operator in eq. (|^) now reduces to: 



pB 



B 



(53) 



In view of the arguments presented in sec. |2.1| above, we find this identification of the two 
distribution amplitudes rather surprising. Indeed the authors of ref. |jT5[ argue that the 
two DAs are likely to have a different end-point behaviour. It would be surprising to us if 
the pseudoscalar matrix element in eq. (^) and the axial one in eq. (TU) could be described 
in the heavy-quark limit by the same distribution amplitude. In this section we examine 
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the weaker proposition in ref. 0; that although 0^ and (p^ are different, they have the 
same contribution to the convolution (0), so that that their difference can be neglected 
when studying B ^ ir form factors. 

In order to assess this claim, we need the B tt form factors ([15|) and (|l^) written as 
convolutions over both longitudinal and transverse momenta: 



and 



^0, 



where 



1 



f 1 



dx / d k\ i d l\T{ 



xrjm? + A;^ 



1 l+^R A(/,,/_l) 



m 



m 



1 - 



7i 



I 



1 



xl+rjm + (A;_L + 
f 1 



+ 



xrjim? + A;^ 
1 



-?7— \& 
m 



B 



(A^x + /I) ■ /j 



(54) 



(55) 



(56) 



with A{ujJa_) = /o"rf£[^^(£,/±) - ^^i£,l±)]. are functions of /+ and and of x 
and A;_L. In the formula for Xq, the contribution in the second line comes from the upper 
diagram of Fig. |^, and the third line from the lower diagram. The second term in the 
formula for X+ comes from the lower diagram. 

The authors of ref. suggest that the value of Fo,+ should remain stable if we replace 
v]/;^ by (or vice- versa) in eqs. ( |55| ) and (|5^). To support this idea, they introduce the 
following "reasonable parametrization" for the distribution amplitudes: 



0?(/+)=5(/+-A)-A5'(/+-A) 



0?(/+)=5(/+-A) 



(57) 



where A is a parameter of order Aqcd- After integrating Xq and X+ over the transverse 
momentum l±, the hard-scattering kernels behave like ln(/+/m) for small /+. Convoluting 
them with (p^ and 0^, we obtain two equal contributions proportional to ln(m/A), but 
with different constant terms. It is therefore argued that the two B meson DAs can be set 
equal to each other in 5 — vr form factors with reasonable precision. 
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The above argument however, raises a number of important questions. First of all we 
note that the two 5- meson DAs in eq. (^7D do not satisfy the equation of motion eq. (Q). 
We also note that this argument implicitly assumes that the dependence on the transverse 
momentum of the modified DAs \E'^ is approximately flat. Of course we accept that the 
model in eq. ( |57D was introduced only for purposes of illustration, nevertheless we do not 
consider it to be sufficiently general to support the case that we can identify (p^ and (p^, 
and use a single DA in the analysis. To underline this point we introduce an alternative 
(and equally unrealistic) model to that in eq. (|57|), but one which satisfies the equation of 
motion ([T^) : 

= - A) and = |e(/+) , (58) 

where A = O(Aqcd), and we define = 1 when < < A, and zero outside 

this range. We now perform the following simple exercise. We take the first term in 
eq. (^) (proportional to \1/:^) and assume (as above) that its dependence on the transverse 
momentum is approximately fiat over [0,A]. We then study the contributions to Xq with 
(f)^ given by each of the two functions in eq. (^5]). In this way we investigate whether it 
matters which of (p^ or (p^ is used in the calculation (at least for this one term). 
If we use (p^{l^) = 0(/_|_)/A, we can perform the integration over and l±: 



f-nfBm —- — / dx d k_i-^^— ^c-{x,k±) , (59) 

l\c Jo J xm? + ki 



where c_ is a function of x and K = k±/A: 
1 



L\l-K'\-S- 2Llog[l + K' + \1- K'\] (60) 



2A2x 

+x - 2x log[2x] + 2x log[x + S + L{1- K"^)] + 2L \og[x + S + L{1 + K"^)] 



and L = A/m and S = ^[x + L'^{1 + K'^)]'^ - AK'^L'^. 

If we use (p^{l^) = 6{l^ — A), we obtain a similar expression, where c_ is replaced by: 

c+ = ^ [log[x + S + L{1- K^)] - \og[2x]] . (61) 

From eq. ( |59D we see that c_ and 0+ have to be convoluted with \I'^/(xm^ + kjj in order 
to compute the contribution of to the form factor Fq{0). When the integration over x 

and k± is performed, the largest contribution comes from the region of small x and k±. 
For small L = A/m, the behaviour of c_ and C+ is given by: 

c+ ~ —IT— (62) 



rn^Lx rn^Lx 



which will be convoluted with the pion DA in eq. (|59D. We see therefore that the model in 
eq. ( |57D was artificially designed in order to yield the same leading contribution ln(AQCD/''^)- 
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Figure 5: Ratio R = , as a function of x, for different values of K = k±/A [A/m = 0.25/5.28]. 



In general, the contributions from (j)^ and 0_ have different leading logarithmic behaviour 
and are significantly different in size. In the model defined by eq. ([58|), the expressions in 
eq. (|62D suggest that (j)^ yields a much larger contribution than c+ and we now check this 
numerically. 

We introduce the ratio R = (c„ — c+)/(c_ + c+) which depends on x and K. The 
identification cj)^ = (p^ = (p^ is justified if R is small in the region which gives the dominant 
contribution to the convolution eq. (|59[) . In fig. ^ we plot the ratio i? as a function of x for 
different values of K = k±/A. For purposes of illustration we take L = A/m = 0.25/5.28. 
We see that R can be as large as 40 % for kj_ < A, even for very small x. The identification 
of the two 5-meson DAs would only be valid for x < 0.05, or kj_ ^ A, i.e. in regions which 
do not contribute very much to the convolution eq. (^). 

The analysis of a different, and perhaps slightly more realistic model, introduced in 
Appendix [A.2| , confirms the possibility of a large error. Of course, we have not estimated 
the actual error on Fq^^ when \E'^ and \E'^ are set equal to each other. In this presentation 
we have considered only a part of the contribution to the form factor -Fo(O), the scale of 
the coupling constant has not been defined, and Sudakov effects have not been included. 
Nevertheless, lessons can be learned from this brief analysis. The hard-scattering kernel 
enhances the small /+ region, where the 5-meson DAs are concentrated. If (p^ and (p^ 
have different end-point behaviours, the difference may be strongly enhanced by the hard- 
scattering kernel. There is therefore no convincing motivation for setting the two DAs 
equal to each other and such a procedure is likely to lead to unreliable results. 

A related issue is the contribution from A{lj) = Jq di {4>^{i) — 0+(^) )• The results of 
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Figure 6: 0{as) corrections to the distribution amplitude for the i?- meson. 



the standard approach, eqs. (|T^ and (|TB|), can be obtained by neglecting the transverse 
momenta in eqs. ( ^Sf ) and (0). We see therefore that a term proportional to A is contained 
in eqs. (|T5[) and (|16|) f\ The contribution of A to the standard formulae eqs. (|T5]) and ([16]) 
is divergent when Z+ — > or x ^ and a large contribution from this region is therefore 
also expected in the pQCD approach. This contribution is dropped however, if one sets 
^+ = (and hence A = 0). 

We conclude that the model approximation in which the two distribution amplitudes 
are set equal to each other is not a generic one and in general will lead to wrong results 
for physical quantities. 



5.2 Derivation of the Sudakov factor for the 5-meson 

In spite of the strong reservations expressed in the previous subsection, let us suppose 
that we can consider a single B-meson distribution amplitude. In this subsection we then 



examine the derivation of the Sudakov factor for the DA of the 5-meson given in ref. |]T0 
and used in subsequent phenomenological analyses. According to the authors of ref. |]lOl , 
there can be an overlap of soft and collinear divergences in gluon exchanges, when /+ is 
much larger than Aqcd? i-e. of order Mb- These divergences can be resummed into a 
Sudakov factor, which is therefore "half" of the Sudakov factor for the pion, since there 
is one light valence quark instead of two. This factor would then lead to a Sudakov 
suppression of configurations with large transverse separations. 

We work in the rest frame of the 5-meson in which the light valence quark has com- 
ponents of momentum of the same order in all directions [O(Aqcd)]- The term "collinear" 
therefore is not meaningful at this stage and no double logarithms are expected to arise, 
as noted by the authors of ref. If one of the light-quark components, /+ say, is much 
larger than Aqcd, the possible Sudakov suppression happens for configurations that already 
highly suppressed. The impact of Sudakov effects is therefore likely to be rather weak. On 
the other hand, in the most likely configuration where all the light-quark components are 
small [O(Aqcd)]; no Sudakov suppression is expected. There is therefore no perturbative 
QCD mechanism to suppress configurations with a large transverse separation but a small 



*This contribution is however not exphcit in eqs. ( JTsI ) and ( |l6| ) because we have used the equation of 
motion eq. to write A(/+) = 
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longitudinal momentum. The consequence of this fact will be investigated in section ^. 

The argument in ref. [0 attempts to proceed along the same lines as were outlined for 
the pion in section H?^, by studying the dependence of the 5-meson DA on /+. We now 



briefly note, without comment, some of the relevant features of this argument, stressing in 
particular the differences with the corresponding derivation for the pion. We then present 
our criticisms of this derivation. 



The authors of ref. ||T0| aim to identify the invariant parameters on which the distri- 
bution amplitude \1'^ = \1':^ = can depend. Because of the presence of the heavy 
quark, many such invariant parameters are available: p"^, p ■ I, p ■ n, I ■ n and n^. It 
is still true that the structure of the gluon propagator in the axial gauge leads to a 
scale invariance of with respect to n, but this argument is not sufficient to deduce 
that depends on a single parameter. 

It is then argued that the eikonal approximation can be used to overcome the diffi- 
culties caused by the presence of the heavy quark. Consider a gluon of momentum q 
attached to the heavy quark line of momentum p — I. When q is soft, or coUinear to 
/, we can use the approximation: 

-7„ = h R . (63) 



{p — I + qY ~ P ■ Q' 

R contains terms that are suppressed by powers of 1/m or that vanish once contracted 
with (J) + m). 

The eikonal approximation leads to a scale invariance of '^^ with respect to p. p 
would only appear in ratios such as (n which are disregarded because 

they cannot lead to a "large scale". The only "large scale" would be /, and the 
relevant invariants would be (/ • n)^ jr? and (/ ■ pY /p^- 

The 0{as) corrections are examined in the particular axial gauge n = (1, —1, 0_|_), in 
which (Fig. |^) is found to exhibit a scale invariance with respect to / in the same 
eikonal approximation. The authors of ref. conclude that (/ • pY /p^ has to be 
discarded and depends only on the scale v' = [l ■ nY/n"^. 

Since \Ef^ depends on a single scale z/', it is possible to replace l+d/dl+ by d/dn°', 
and to derive a differential equation for the 5-meson distribution amplitude: 

/+^v&^(/+, 6; /.) = ^ [/C(6/., (?(/i)) + g{u V/i)] vI/^(/+, 6; /i) . (64) 



This equation can be solved using the same techniques used for the pion in section 
leading to: 



^^(/+,6;/i) = e-^'^^'^'^^'^exp 



0«(/+;/i = l/6)+O(a,(l/6)), (65) 



where s is defined in eq. (|39|). We can rewrite the Sudakov factor in terms of Z_| 
s(l, h,v') = s{l+/m, b,m/r), where s is defined in eq. 
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We now express our concerns about this derivation, which include the identification of 
the two 5-meson distribution amphtudes, the use of an eikonal approximation for the heavy 
quark and the determination of the scales on which the 5-meson distribution amplitude 
depends. 

• It is not clear to us that multiplicative, RG-invariant Sudakov factors arise separately 
for \E'-^ = {"^f + \E':^)/2 and = {"^f - \E':^)/2 (the latter is neglected in ref. p^). 



• The eikonal approximation eq. (|6^) is relevant when /_|_ is soft. But the resummation 
of double logarithms is performed in the region where /+ is large. The relevance of 
eq. (|63|) in this case is not clear to us. 

• A demonstration of the /-scale invariance of the 0{as) corrections to in a partic- 
ular gauge, does not ensure the scale invariance of the function \E'^. And even if this 
property was satisfied by \E'^, {I ■ p)"^ /p^ and (/ ■ n)^ jr? are both /-scale dependent 
and therefore could not be relevant variables for . If \E'^ is scale independent with 
respect to n, p and /, it cannot depend on any "large scale" . The invariants with this 
property are 

(p-n)V[pV] = £(l + ^y , [l-pfn^l\p\l-nf\ = y (66) 

where n = {\/^/2, p/y/2,0±) and only the /+ component of / has been retained. 

Moreover, the questions raised by the derivation of the pion Sudakov factor remain to 
be answered: 

• How is the gauge dependence of \E'^ cancelled by that in the hard-scattering kernel ? 

• Is it legitimate and sufficiently precise to identify (pB and \E'^ ? 

We have tried and failed to understand the derivation of the Sudakov factor for the 
distribution amplitudes of the 5-meson. In our view it is based on a number of unjustified 
assumptions and steps. Given that the dominant contributions come from the region of 
small /+, it seems to us that Sudakov effects are unlikely to be significant for the B- 
meson. If this is the case however, the integrals in eqs. (|l^) and ([T6|) may give significant 
contributions (or even diverge) at small /+. 



6 Numerical study 

In the previous sections we have explained our reservations about the pQCD formalism 
used to evaluate the semileptonic form factors and other physical quantities. We now tem- 
porarily set aside these reservations and perform a numerical study of the pQCD evaluation 
of the 5 — > TT semileptonic form factors. We are particularly interested in two questions: 
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1. Are the results sensitive to the choice of distribution amphtudes? 

We will see that the answer to the first question is yes, so that the form factors are 
not calculable in practice. 

2. For "reasonable" choices of distribution amplitudes, does all of the contribution come 
from the perturbative region? 

The answer to this question is in general no, so that the form factors are not calculable 
also in principle. This conclusion is in addition to our theoretical reservations. 



6.1 pQCD expression of 5 ^ tt form factors 



The formulae for the form factors in eqs. 
parameter space: 



and (|56D can be rewritten in impact 



dx I dl+ / b±db± 



J c±dc_i 



^7r(a;, b±; fj.)Ho^+{x, /+, b±, c±, r], m; c±; fi) 



where 



,7lCl. 
1^ 



'1 + xri)G(xl+rim, b±, xrjm , c±) H Gil+rjm, c_l, xl+rjm, 6j 

m 



and 



fnfBm- 



lv7 



[1 + xri)G{xl-^rim, b±, xrjm , c±) G{l+rim, c±, xl+rjm, b±) 



m 



with the function G given in terms of Bessel functions by: 

G(A&±,5,cx) = Ko{^b^)[e{b^-c^)K^{^b^)h{^c^] 



+e{c^ - b^)KoiVBc^)Io{VBb^)] 



(67) 
(68) 



(69) 



(70) 



The angular integrations in the transverse plane have been performed, since the distribution 
amplitudes are independent of the orientation of the transverse momenta. 

There are higher order corrections to the hard-scattering amplitude of order log(t//i). 



where t is the large scale associated with the hard gluon. In ref. [10| this was taken as 



t = max(v^rj~CT7m, 1/c_l). In order to avoid large corrections, we can perform the 
evolution of the hard-scattering kernels ifo,+ from fi to t using the RG equation: 



^ d d\ 



(71) 
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where 7^ = —asj'R is the anomalous dimension of the quark's wave function in the axial 
gauge. For H^ j^if) we take its expression at lowest order in a^. At this order H does 
not contain a gauge dependence able to cancel the gauge dependence of and . In 
particular, there is no resummation of soft gluon contributions in the hard-scattering kernel 
in a similar way to the work of ref. for hard elastic scattering of mesons. We obtain: 



fnfBm 



,7lCf 



dx / dL b\ db 



J c±dcj 



xrj 



N, Jo 

x0^(x; 1/c±) exp[-S'(x, l+,r], m, b_i_, c±)]as{t) 

'1 + xri)G(xl+rim, b±, xrirn? , c_l) + -^GU^rjui, c±, xl+rjm, bj_) 

m 



(72) 



and 



dx I dL bi db 



J c±dc_i 



X 



x0^(x; l/b±)(j)^{l+; 1/c_l) exp[-S'(2;, /+, 77, m, b±, c±)]as{t) 

[1 + xr])G{xl+r]m, b±, x'qm^, cj_) —G{l^r]m, c±, xl+r]m, 6j 

m 



(73) 



where the function S combines Sudakov factors and evolution-related terms: 
S 



s{x, 6_L, rimr) + s{x, b±,rimr) + s{l^/m, cj_, mr) 



(74) 



+2f %,igm + ^[] -iMn)). 



The function s is defined in eq. (@). We follow refs. [[l^, ^ and take the particular values 
Gi = 1, C2 = \/2 and r = 1, so that 



g = log[MQ/(y2AQCD)], and 6 = log[l/(6AQCD)] • 



(75) 



In eqs. (|7^ and (|73D, the transverse separations are integrated from to 1/Aqcd, 
because the Sudakov factors suppress any configuration with larger transverse separations. 
We can accept this statement in the case of the pion (6_l). But the Sudakov factor for 
the i?-meson is weak for the usual configurations where Z+ is of order Aqcd or less. We 
may therefore worry that the integration with respect to cj_ is artificially cut-off, and omits 
regions which might yield significant contributions to the form factors (such as small /+ 
but large c_l). 



6.2 Models of distribution amplitudes 

In order to compute the B —>■ n form factors, we have to choose models for the distribution 
amplitudes of the pion and the B meson. It can be seen from eq. (^) that the modified 
pion DA '^T^{x,b; fi) is the product of three factors: a Sudakov exponential, an evolution- 
related term and the BLER pion DA 0^(x;/i = l/b). In ref. |T^, the 6-dependence due 
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Figure 7: Models (f>fi), 4>f2) ^^^^ for the B distribution amplitude, as functions of ^ = ^+/m. 



to the renormalization scale of (p-^^x] jj, = 1/b) has been neglected in eqs. ([7^) and (|73|), 
because its numerical impact is believed to be small. We will follow this prescription here. 
In a similar way, the 6-dependence of the S-meson DA in eqs. (|72D and ([73| ) is neglected. 
We use two different models for the leading-twist pion DA [0, : 

(Pi^'^ (x) = Qx{l~ x) [1 + aaC^^/^^ (2x - 1) + a4Cf {2x - 1)] , (76) 
= 30x(l-x)(l-2x)2 , (77) 

where C^^'^\u) are the Gegenbauer polynomials. Even though we neglect evolution effects 
in the computation of the B ^ n form factors, the BLER distribution amplitudes do 
depend on the renormalization scale fx. In model (p^^ the /i-dependence is contained in the 
coefficients a2 and a^, model (j)^^ was proposed originally at the scale /zq =0.5 GeV [O . 



We will use three models for the "common" B-meson DA; in all three models the 
contribution from the region with /+ > m is suppressed. Defining ^ = /+/m, the first two 
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models are 



proi 



'/'fl)(^H 
</'f2)(^H 



iV(i)V^(l-Oexp 



m 



2u' 



and we introduce a third model, with a different end-point behaviour: 



0f3)(^H 



iV(3)(l 



(78) 
(79) 

(80) 



The normalization constants A'"(j) are obtained from the integral / dl^(j)^{l^) = 1. Fig. 
shows the general shape of the three models. /+ is concentrated around uj/\^ in the first 
model and we expect a; = O(Aqcd)- The second model has a broad distribution in /+, and 
may therefore be considered highly unlikely to be physical. The third model is concentrated 
in the small /+ region. 

We stress that the three models in eqs. (|78D-(pOD are introduced to study the dependence 
on the expressions (^) and (|73| ) on the end-point behaviour of the distribution amplitudes. 
Our reservations about the consistency of using a single wave function for the i?-meson 
remain of course. In particular the question of whether 1/ is finite for each of the three 
models depends on whether the functions in eqs. ([78|)-(|80|) are interpreted as corresponding 



to 0^ or 0_, and cannot be answered in the approximation in which a single distribution 
amplitude is used for the S-meson. 



6.3 Dependence on the shape 

We now use eqs. ( [72D and (|73D to compute the form factors and Fq for the models of 
the 71 and 5-meson DAs introduced in section We are able to reproduce the quoted 
results in ref. for the first two models for the B meson 0^-) and 0^-,. We wish to 



study the model-dependence of these results. For illustration we will present our results 
for the favourable kinematic situation rj = 1 (i.e. = 0, where Fq = F+). The results 
are presented for the following choice of parameters: fs = 0.19 GeV, m = 5.28 GeV, and 
Aqcd = 0.25 GeV. 

For the 5-meson DA we consider the three models introduced above and vary the shape 
parameter u. We use the two models for the pion DA in eqs. (|76D and (^). Eq. ([7^ ) is 
linear in 0^. If we choose the expansion in Gegenbauer polynomials (f)^^ = (p^^ we can 
give separately the numerical contributions from each of the three terms in eq. ([7BD, i.e. 
the contribution from the asymptotic DA and those proportional to a2 and 04. The corre- 
sponding values are indicated in Tab. |^. For instance, if we choose (j)^ = 0^^, u;=0.4 GeV 

and 0^ = (f)i''\ the form factors at = are Fo,+ (0) = 0.12 + 0^2 ■ 0.25 + 04 ■ 0.25. 

We see that the computed values of Fo,+ (0) depend very significantly on the distribution 
amplitudes of both mesons. Models concentrated around small values of /+ (0^-) for small u) 
or X (0^^-*) give larger values than models which have a broader spread (0^), asymptotic pion 
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5-meson DA 


Pion DA 


iViO QfJl 


UJ 


Ac 






P 7 




n 1 

U. 1 


n 9S 


n zLfi 


n Af, 


n 

u.oo 




n 9 


n iQ 


u.oo 


u.oo 


n 49 




0.3 


0.14 


0.29 


0.29 


0.33 


n A 


n 1 9 


n 9^; 


n 9^; 


n 98 
u.zo 




U.O 


n 1 n 


n 99 


n 9fi 


n 9^; 






n DQ 


n 1 Q 


n 9'? 

u. z,o 


n 99 




97 ^ 




n 1 9 




U. LO 




-25.5 


0.04 


0.10 


0.13 


0.10 




-22.5 


0.03 


0.09 


0.11 


0.09 


-20.5 


0.03 


0.08 


0.11 


0.09 




-17.5 


0.03 


0.08 


0.10 


0.08 




0.1 


0.36 


0.52 


0.51 


0.70 




0.3 


0.21 


0.36 


0.38 


0.45 


<^f3) 


0.5 


0.16 


0.29 


0.32 


0.35 




0.7 


0.14 


0.25 


0.28 


0.30 




0.9 


0.12 


0.23 


0.26 


0.27 



Table 1: Dependence of Fo{0) = F+(0) on the B -meson and pion DAs. The models of _B-meson DA 
depend on the shape parameter uj (in GeV for models 1 and 3, GeV^ for model 2). The columns "As.", 
"q!2" and "0:4" correspond to = 6x{1 — x), (f)^^ = 6a::(l — a::)C2^^^'' (2a::— 1) and (/)^ = 6a:;(l — a;)C4^^^^(2a;— 1) 
respectively, ^t"' is a linear combination of these three terms. "C-Z" corresponds to the model 

DA). This can be readily understood since the hard-scattering kernel enhances distribution 
end-points (without invoking Sudakov effects the contributions from the end-point regions 
are divergent). 

The conclusion of our investigation is that very good control of the behaviour of the 
DA's at the end-points is needed in order to be able to compute the form factors with a 
precision which would be useful for phenomenological studies. From table |l] we see that 
with our present knowledge of the distribution amplitudes this is not the case. 

6.4 Dependence on the cut-ofF in impact parameter space 

In the previous subsection we have seen that with our present knowledge of the mesons' 
distribution amplitudes, it is not possible to calculate the form factors with the required 
precision. We now investigate whether, for model distribution amplitudes of the form 
used above, all (or almost all) of the contribution to the form factors comes from the 
perturbative region of phase space. In the preceding sections we have seen that in order 
to set "^nix, b±, Q; fi = l/b±) equal to (p-^ix; fj, = l/b±) (see eq. (^) ) we require b± to be 
small p3 (there is a similar requirement for the 5-meson). We now check whether all of 
the contribution does indeed come from the region of small impact parameters. 
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Figure 8: Dependence of Fo,+ (0) on the cut-offs and for the integration over b± (pion impact 
parameter) and c± (B- meson impact parameter). For purposes of iUustration, in this figure we choose 
as the model distribution amphtude for the B-meson, and the asymptotic distribution amphtude for the 
pion. 



In fig. |g we evaluate the integrals in eqs. ([7^ ) and (^), but with a cut-off introduced 
for the impact parameters for the pion or the i?-meson. For purposes of illustration we 
take for the 5-meson's DA, and the asymptotic distribution amplitude for the pion, 
but similar plots can readily be obtained for the other models of distribution amplitudes. 
The figure contains plots with the integrals over the impact parameters performed over 
the regions < b± < and < cj_ < 1/Aqcd (left-hand figure) or < 6j_ < 1/Aqcd and 
< < 6^ (right-hand figure). Fig. ^ shows the dependence of Fo^+(0) on the cut-offs 

(left) and b^ (right). In order for the calculations to be consistent, we require that the 
contribution from the regions of large impact parameters is negligible, and hence that the 
curves in fig. ^ reach a plateau for values of 6^ and 6^ in the perturbative region. 

The purpose of our investigation is to check whether this is the case. In general, and for 
the models used in fig. |^ in particular, the answer is clearly no. Even if we optimistically 
take 500 MeV as the value above which perturbation theory holds, we see that the curves 
in the figure are far from saturating at this scale. There is a significant (and uncalculable) 
contribution from the nonperturbative region of phase space, where eqs ( p5D and ( ^2]) 
are affected by large corrections, impossible to estimate. Contributions from regions with 
larger impact parameters cannot be calculated reliably and hence we conclude that pQCD 
calculations for the form factors are not valid. 



A similar study was reported in ref. [T0|, in which the authors impose that for a con- 
sistent pQCD computation, most of the contribution should come from the region where 
as{l/b^)/TT and as{l/c±)/7r are smaller than 0.5, i.e. the impact parameters are smaller 
than b^ = b^ = b'^ < 0.6/Aqcd- The authors concluded that pQCD approach was rela- 
tively self-consistent, since a large contribution comes from the perturbative region. Our 



discussion above makes it clear that we do not accept the conclusion of ref. ^0 



The authors of ref. [0 considered a weaker criterion, that most of the contribution 
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Figure 9: Contribution to Fo_+(0) as a function of the coupling constant as{t) in the hard-scattering 
kernel. For purposes of illustration, in this figure we choose i/)^^ as the model distribution amplitude for 
the B- meson (with lo — 0.1 GeV), and the asymptotic distribution amplitude for the pion. 

should come from the region of phase space in which the maximal available virtuality (lon- 
gitudinal or transverse) of the gluon, t = max{y/ 1 +xrim, l/b±, l/c±) is in the perturbative 
regime. They take t as the scale of the coupling constant in the hard-scattering kernel. 
We have argued above that this condition is insufficient for the consistency of the cal- 
culations, l/b± and l/c± must also separately be in the perturbative regime. Of course, 
since t > l/b± and t > l/c±, for any set of distribution amplitudes a larger fraction of 
the form-factor comes from the region is which t > where is a cut-off defining the 
perturbative region of phase space, than from the regions in which l/h±_ > or 1/c_l > f^. 

We have however considered this weaker criterion, and studied the contribution to the 
form factors coming from the region of phase space where ai < as{t)/7T < a2- Fig. ^ is 
the resulting histogram for (p^ = 0^^, uj=0.1 GeV and (p^^ = The last bar on the 

right-hand side is the contribution of as{t)/7i larger than 0.9. We find that the fraction 
coming from the region in which t is not perturbative, is large. Numerical studies of other 
models for the distribution amplitudes indicate that the contribution from the region of 
non-perturbative t can be large, but this depends very much on the choice of distribution 
amplitudes. It may be relatively small for some choices of the distribution amplitudes [0 
and large for others. 

Sudakov suppression of configurations with large transverse separations is not efficient 
enough for i? ^ tt form factors, whereas it worked reasonably well in the case of the pion 
electromagnetic form factor [^. This should not surprise us. Recall that there is a 
fundamental difference between the two processes in the BLER approach: the first suffers 
from long-distance divergences, whereas the soft contributions are finite in the latter case. 
One would therefore expect that a much stronger Sudakov suppression is required, and not 
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achieved, for a self-consistent pQCD approach of S-decays. 

As explained above, we have serious concerns about the derivation and the expression 
of the Sudakov factor for the 5-meson. Note that the weakness of the Sudakov suppression 
in eq. (0) can be expected on general grounds for the i?-system. We have seen that a 
Sudakov effect was argued for i?-mesons with a large longitudinal light-quark momentum 
/+, which is a highly unlikely configuration. On the other hand, there should be no effect for 
the standard situation: /+ ~ /_ ~ ~ Aqcd- "Large" transverse separations {c± > Aqcd) 
are not suppressed by this mechanism. Therefore, even if an expression of the form ( ^2]) 
could be derived for the Sudakov factor in the i?-meson, it would lead to large contributions 
coming from nonperturbative regions. 

We conclude this section by restating that the pQCD predictions for B ^ it form 
factors receive a substantial (but uncalculable) contribution from configurations with large 
transverse separations. The Sudakov suppression appears to be too weak to cure this 
problem. Although the numerical details of our study depend on the choice of models for 
the distribution amplitudes, the qualitative conclusion is general. We are therefore forced 
to conclude that the pQCD predictions for B ^ n semileptonic form factors (and other 
physical processes which are affected by end-point singularities) are not valid. 



7 Conclusion 

In this paper we have studied Sudakov effects in predictions for B —>■ tt semileptonic form 
factors. In the standard approach, these form-factors are not calculable due to the pres- 
ence of long-distance effects from end-point regions of (longitudinal) phase-space. We have 
investigated the claim that Sudakov effects suppress these long-distance contributions suf- 
ficiently for the form-factors to be calculable reliably and precisely in perturbation theory. 
Our conclusion is that this is not the case. The same arguments can be used to conclude 
that Sudakov effects cannot be invoked to make reliable predictions for other processes 
which have end-point singularities, such as power corrections to the amplitudes for exclu- 
sive two-body 5-decays p, ||, |^, |31]] . 



Among the reasons for our conclusion are: 

1. As explained in detail in section ^ it is not possible for us to accept that the formalism 
currently used to derive the Sudakov factor for the 5-meson is theoretically sound. 

2. Even if one accepts the Sudakov factors for both B and vr mesons, the uncertainty in 
the mesons' distribution amplitudes (particularly at small x and /+) means that the 
form factors for B irlui decays cannot be evaluated with sufficient precision to be 



phenomenologically useful. This was investigated numerically in section 3.3 



For these i?-decays the kinematic parameters are such that substantial contributions 
to the form-factors come from the nonperturbative region (see sec. |6.4| ), and are 
hence uncalculable. Sudakov effects are too weak to suppress contributions from the 
regions of phase space with large transverse separations. We therefore conclude that 
the pQCD approach is invalid for semileptonic -B-decays. 
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Following the completion of this work, there appeared the interesting paper by H. Kawa- 



mura et al. , in which the heavy quark theory and equations of motion are used to show 
that, under the assumption in which the three-parton distribution amplitudes are neglected, 
the two B-meson distribution amplitudes (p^ and (p^ can be determined (more generally 
it is shown in this paper that the two- and three-parton distribution amplitudes can be 
related). These results reinforce our conclusions, in particular it is shown in ref. [0 that 
(f)^ does not vanish at the end-point. Numerical studies with these distribution amplitudes 
lead to the same conclusions as presented in section |^. 

As mentioned in the introduction, our inability to evaluate the power corrections in two- 
body nonleptonic -B-decay amplitudes {B M1M2, for two mesons Mi and M2) severely 
limits the precision with which we can deduce fundamental information from experimental 
measurements of the rates and CP-asymmetries. Our conclusion, that Sudakov suppression 
of long-distance effects is too weak (and too unreliable) to be useful in extending the 
range of applicability of perturbative QCD in 5-physics, is therefore a disappointing one. 
Nevertheless, given the fact that many phenomenological studies of 5-decays are being 
performed which are based on the Sudakov suppression of end-point singularities, we felt 
that it was important to articulate our concerns about the validity and reliability of this 
approach. 
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A A separable model for the B-meson 

In this appendix we introduce a model for the two leading-twist S-meson distribution 
amplitudes (in longitudinal and transverse momentum space). This model is introduced 
for illustrative purposes and is used in section |^. It satisfies the constraints arising from 
the equation of motion. 



A.l Constraints from the equations of motion. 

We will generalise eq. (^) to include the effects of transverse momenta. Since the hard- 
scattering kernel is independent of /_, we can consider eq. (|]) for = (but 0). The 
most general decomposition is: 



[0Mz)bo,mB{p)) 



B 



^ + 171 



2^'l{z',t) + 



^^{z^,t)-^f{z^,t) 
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with m = Mb = mb, p = mv and t = v ■ z {v = p/Ms)- A path-ordered exponential is 
imphcitly present in the gauge-independent matrix element. We can introduce the Fourier 
transforms: 

(P'Zx_ dz. 



(27r)2 27r 



.g«{'+^--'±-^±)v[/ 



B 



82 



We have defined / = (/_|_/v^, 0, and z = (0, Z-\/2, z±), so that z"^ = —zl and t = z_. 

In ref. |16| constraints were derived on the leading-twist 5- meson DAs, neglecting the 
effects of three-particle {qqg) and higher Fock states: 



9§2 



t 







22=0 



2,T,B 







(83) 

dz^ ' A dt^ " ' ^^^^ 

where depend on z"^ and t = v ■ z. We consider a model in which the dependence on 
the longitudinal and transverse momenta is factorized: 

*f(^+,rx) = 0|(Mxrf(/x), (85) 

where we chose the normalization conditions 



dl+^l{l+) = 1 and / d^l±T±{l^) = 1 . 



In this case, the two constraints become in the momentum space: 

B '^'^^ 



dL 



(86) 



^7) 



where = / d'^l±_?^T+{lj 
tion: 



Combining the two equations, we obtain the differential equa- 



and the corresponding (normalized) solutions: 



2A2 



and 



2 / 



2A2 



^9) 



(90) 



The value of A depends on the model used for the transverse momentum, and measures 
the dispersion of its values. For example, the step-function distribution: 

_ f l/(7rAi) for 0<1^<X^ 
I otherwise 



r+(/j 



(91) 



leads to A = \±/V^, whereas the Gaussian distribution: 



27rai 



exp 



2a 



(92) 



has A = V^a_L. 
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A. 2 Can we set = ^^^7 

As an exploratory exercise we estimate the error in the calculated value of the form-factors 
(see eqs. (H), ( ^5|) and (|^)) caused by setting \E'^ equal to For this exercise, we 
take for the DA of the i?-meson the separable model of App. [A.l| , with a fiat transverse 
distribution, eqs. (^) and (pT]). We set A_l = A, with A = O(Aqcd)- For this model the 
first term in eq. ( ^5]) is: 

pi poo /• /• 

S_{7]) = / dxMx) / dl+^^il+) / (fk^T^k^) / d^LT+{l^) (93) 

1 1 

XT]- 



xl^rjin + (A;_|_ + l^Y xrjrn^ + k\ 

27]- \ 2 ^1 ^oo /-A /-A 



/ dxM^) rf/+0^(/+) / A / (94) 
\ttA^/ Jo Jo Jo Jo 



1 



X7]m^ + kl ^{xl+r]m + kl + llY - Aklll 

r, rl roc rk^ rlV^ 

/ dx(f)^{x) / dl+<p^{l+) / dn dX (95) 
JO Jo Jo Jo 



1 1 



xvm'^ + 1^ ^{xl^7]m + K + Xy -AkX ' 
The following integral is useful for the angular integration over the transverse momenta: 

/ -dd = , for a > 6 . (96) 

Jo a + bcosO - 

In eq. (0), the integrals over A, k and /+ can be performed numerically and we find : 

I.ix,v) = Ti / rf/+0^(/+) dK dX (97) 
A* JO JO JO 

1 1 

^Vm^ + 1^ ^{xl+T]m + K + A)2 - 4kX ' 

The same integral can be computed with \E':^ instead of \E':^, and the result is called /+(x, rj). 



For illustration, the two functions are plotted in Fig. [TO for the case i] = 1. Throughout 
the region in x, J_ remains about twice as large as /+. 

In order to estimate the error on the form factors we need to take a model for the DA 
of the pion. For example, if we take the asymptotic form, 07r(x) = 6x(l — x), we obtain 
S_{r] = 1) = 1.32 GeV"^ and S+{r] = 1) = 0.59 GeV-^ The identification of ^+ and ^_ 
leads therefore to an error of at least 30% in this comtribution to the form factors. 
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Figure 10: First term I-{x,rj = 1) in eq. (55), with the separable model for '5;^, and corresponding 
integral I+{x, rj = 1) when is replaced by "if^. 
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